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Abstract 

In this paper we study a problem in radiotherapy treatment planning, 
where the evolution of the radiation field is governed by a deterministic 
Boltzmann transport equation. We show existence, uniqueness and reg- 
ularity of solutions to an optimal dose distribution problem constrained 
by the Boltzmann Continuous Slowing-Down equation in an appropriate 
function space. The main new difficulty is the treatment of the stopping 
power term. Furthermore, we characterize optimal controls for problems 
governed by this transport equation. 

1 Introduction 

Besides surgery and chemotherapy, the use of ionizing radiation is one of the 
main tools in the therapy of cancer today [31 HI El E]. According to WHO data, in 
the year 2007 there were about 11.3 million new cancer cases. More than half of 
the patients that are treated receive radiation therapy at one point during their 
treatment. Since the early days of radiation treatment high energy photons have 
been the most important type of radiation. Other types of radiation include high 
energy electrons and heavy charged particles like protons and ions. The latter 
type of radiation is of growing importance, but has not reached the widespread 
use of photons and electrons, yet. The aim of radiation treatment is to deposit 
enough energy in cancer cells so that they are destroyed. On the other hand, 
healthy tissue around the cancer cells should be harmed as little as possible. 
Furthermore, some regions at risk, like the spinal chord, should receive almost 
no radiation at all. 

It is still current practice that treatment plans involve several fixed beam 
directions which arc selected by an experienced physician by hand. Radiation 
facilities where the beam head rotates around the patient and where the beam 
is shaped by multileaf collimators are entering clinical practice. These methods 
have become known as Intensity- Modulated Radiation Therapy (IMRT). Pa- 
tient motion during treatment is also one of the future challenges in the field 
of external beam radiotherapy. For instance, tumors near the lung move due 
to breathing. Techniques addressing this problem have become known as 4D 
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radiotherapy (4DRT) [3], meaning that time, the fourth dimension, also has to 
be taken into account. A further technique, named Image-Guided Radiotherapy 
(IGRT) is currently being developed. In this method, the radiation is used to 
create patient images during treatment. All of these novel techniques require 
mathematical modeling and optimization techniques [141 123j . 

Before the treatment of the patient can be started, the expected dose dis- 
tribution, i.e. the distribution of absorbed radiative energy in the patient, has 
to be calculated. Most dose calculation algorithms in clinical use rely on the 
Fermi-Eyges theory of radiation. In recent work |16j . however, it has been 
shown that these can produce errors of up to 12% near inhomogeneities. 

This work is based on dose calculation using a Boltzmann transport equa- 
tion [T^l [U] . Similar to Monte Carlo simulations it relies on a rigorous model 
of the physical interactions in human tissue that can in principle be solved ex- 
actly. Monte Carlo simulations are widely used, but it has been argued that a 
grid-based Boltzmann solution should have the same computational complexity 
[nil]. Furthermore, Monte Carlo can only be used in derivative-free methods 
for optimal dose distributions. In constrast, when optimizing using Boltzmann's 
equations it is possible to exploit structural information for numerical and an- 
alytical purposes of the optimization problem [251 Ell ESI ■ 

Our starting point is the Boltzmann equation for particle transport in a 
medium: consider a part of the patient's body which contains the region of the 
cancer cells. We assume that this part of the body can be described as a convex, 
open, bounded domain Z in M^. Furthermore, we assume that Z has a smooth 
boundary with outward normal vector n. The direction, into which the particle 
is moving is given by e S*^, where S*^ is the unit sphere in three dimensions: 



Here, -0 can be thought of as being the number of particles at x £ M'^ with 
energy e, and direction $7 g 5*^. Scattering is determined by the total cross 
section at and by the scattering kernel CTs, which can be seen as the probabilty 
that a particle with initial energy e' and initial direction fi' has energy e and 
direction after the scattering event. 

For high energy particles, small angle and energy changes are very likely, thus 
the scattering kernel s is very forward-peaked. This fact is utilized to derive 
the Boltzmann Continuous Slowing-Down (BCSD) approximation |18| . This 
model still allows large-angle scattering (which is important in radiotherapy 
applications) but describes energy-loss differentially. 




a, (cc, e', e, Vl! ■ n)^jj{x, e', Vt')dVt' de' - at{x, e)'4>{x, e, Vi). 



-Tr{S{x, e)i;{x, e, il)) + Q ■ V,^<a;, e, n) 




as{x, e, Vl' ■ Q,)ip{x, e, Q,')dil' — at{x, e)ipix, e, ft), 

(1.1) 



where 
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and the stopping power is 



S'(x, e) = 27r / / {e — e')(Ts{x,e,e' , ^)dijde' . 



This equation can be viewed as an initial boundary-value problem for 
{x,€,fl) e Z X (0,00) X S^, being the unit sphere in R^. To formulate 
boundary conditions, we define the in- and outgoing boundaries as 

r± := {{x, n) edZ X : n{x) ■ n > (<)0} 

and prescribe 

^(x^n) ^ q{x,n) on r_. 

The "initial condition" is 

ip{x, Qo, rt) = 0, 

meaning that there are no particles with arbitrary large energy. 

Different other approaches exist. For a review on neutral particle codes that 
have been applied to the dose calculation problem we refer the reader to [H] . 

A number of functionals and methods have been devised to describe the 
effect of radiation on biological tissue, cf. the extensive lists of references in the 
reviews [8] and [23]. It is clear that the amount of destroyed cells in a small 
volume, be they cancer or healthy cells, is not directly proportional to the dose 

D{x)= [ [ S{x,e)^{x,e,u})dujde (1.2) 

deposited in that volume. However, no single accepted type of model has 
emerged yet. Moreover, current biological models require input parameters 
which are not known exactly This is why the authors of [13] opted not to 
investigate these models but rather to focus on some general mathematical cost 
functionals. A quadratic objective function together with nonlinear constraints 
was identified as the most versatile model. Divide the domain into tumour tis- 
sue, normal tissue and a region at risk: Z ~ Zt U Z^ U Zji. We prescribe a 
desired dose distribution D, which usually has a constant value in Zt and is zero 
elsewhere. The problem of optimal treatment planning is to find an external 
beam distribution q such that 

J=^[ (D-Drdx + ^l [D-Dfdx + ^j [D-Dfdx 

is minimial. Additionally, we might add the constraints 

D > Dnun in Zt and D < Dmax in Zr, 

which ensure that all of the tumour tissue is sufficiently irradiatied and the 
region at risk receives a limited dose. 

2 Main result 

The target area is modelled by a bounded, convex domain Z cM.^ with smooth 
boundary dZ. All results extend to the case Z C K" and S"~^. 
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We make some additional assumptions in order to provided a concise math- 
ematical treatment. The first assumption is no major restriction. We assume 
that there is a maximal energy denoted by e^aax- Later, we will introduce a 
transformation e' = Cmax — e and solve an initial-boundary-value problem with 
initial values prescribed for e' = 0. 

Our second assumption is that the averaged scattering coefficients at and 
(Js do not depend on energy. This means that the elastic part of the scattering 
process is independent of the energy of the incident particles. This is not satisfied 
for elecron scattering. Our main purpose here is to deal with the stopping power. 
Thus, we postpone the treatment of energy-dependent scattering coefficients to 
future work. The stopping power itself is assumed to be independent of space. 
A possible reasoning for this assumption is that dose calculations are performed 
with data coming from a voxel-based CT scan. The stopping power is constant 
in each voxel and we couple the solutions to the transport equations over the 
different voxels. 

Under these two assumptions we derive rigorous results on existence, unique- 
ness and regularity on optimal controls. We are thus interested in solutions to 
optimal control problems on Z x S*^ x [0, Cmax] subject to the Boltzmann con- 
tinuous slowing-down equation in the following form 

d^S{t)^ + nV^^Pix, e,n)+at {x, n)ij{x, e, n) = (2.3a) 

/ as{x,n' ■n)i:{x,e,n')d^}' + q{x,e,i}), (2.3b) 

ip{x, n,€)^0 on r_, V(a;, n,0) = on Z X S^. (2.3c) 

We denote = max(j;, 0) and we introduce the function spaces 

iL = e L\Z xS^x [0, e^ax]) ■■q>0 a.e. }, 
D{A) = {if e L^{Z X S^) : -nV^ip e L^{Z X S^), ip\T_ =0}, 

and 

W2 = {V e L^Z xS^x [0, e„,ax])); ^ + f7 • G i'(Z X 5^ X [0, e„,ax]); 
V'(-,-,0)eL2(Zx52);V'|r_x[o,e._] ei2(r_ x [0,e,„ax],f^-n(x)dr_)}. 

We impose the following assumptions on the coefficients: 
Al The fimctions at and as are non-negative a.e.. 
A2 The functions satisfy at £ L°°{Z) and cr, e L°°(Z x [-1, 1]) 
A3 The scattering kernel is uniformly bounded for all a; G Z : 




as{x, ^)d^ < c. 



A4 The stopping 5* — 5(e) is a strictly positive and continuous function. 
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electron rest energies 



Figure 1: MoUer stopping power for water as a function of energy. 



Remark 2.1. The assumptions (A1)~(A4-) are satisfied for any reasonable scat- 
tering and absorption kernel. Figure [7] shows the stopping power for Moller in- 
elastic scattering for water. In units of the electron rest energy, it can be written 

as m 



S{e) = 



2^rlp{e + 1)^ 



2 In 



1 



2e + 1 



■In 2 



2(e+l)2^ 4 ^ 

Here, eb is the electron binding energy in a water molecule. 

The positivity and smootheness of S afow us to introduce yet another new 
variable e = r(e) as the unique sohition to 



1 



dr 



r(0) 0. 



Then the quantity ip{x, e, ft) := 'ipi^i ''"^ ^) satisfies a transformed transport 
equation without stopping power term, but with a modified right hand side. By 
investigating this transformed equation, we obtain the following main result: 

Theorem 2.2. Assume (Al)-(A4). Letui £ L°°{Z) be positive a.e., let q ^ L^^ 
and let Tp G L'^{Z x S'^ x [0,emax])- 
Then, the problem to minimize 



ai{x) 



'[0,emax]xZ 

subject to 



4' — ipdfl ) dxde + 

S2 / J[0,e„ex]xZxS2 



a2{x) 



[q — qYdxdVlde 



GsijjdVL' + 5, 



S2 



il){x,n,e) =0 on r_, ^{x,n,Q) = Q on Z xS"^ 
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admits a unique weak solution {ip* , 9*) £ (^"^([O, Cmax], L'^ {Z x S^)) x L^j^. Under 
the assumption ijj e W2 and ip* G W2 H C''([0, emax], ^(^)), 'is a local 

minimum, if there exists a weak solution X* g C'^([0, Cmax]? D{A*)) fl W2 of the 
first-order optimality system: 

d^Sip + nVtp + cTtV' = / (Tsipdn + q, (2.5a) 

^p{x,fl,e) = on r_, 'ilj{x,i},0) ^ on Z x , (2.5b) 

-Sd^X-n\7X + atX^ I asXdn' + ai [ {ij - 4>) dn, (2.5c) 

X{x, n,e)^0 on r+, X{x, n, e^ax) = on Z X S^, (2.5d) 
q* = {q* - X* - a2{q* - q))+ . (2.5e) 

3 Existence and uniqueness of the minimizer 

The proof of Theorem 12.21 is spht among several results. We start with the 
notion of a weak solution. 

Definition 3.1. Let T be a strongly continuous semigroup of operators on a 
Banach space X with infinitesimal generator A. Let f S L^{[0,T], X), T > 0, 
S : R ^ IR+ be a continuous, positive real-valued function and let r be the 
solution to ^ — -gj-^ with r(0) = 0. 

We call ip G C'^{[0,T],X) a weak solution to 

dtip = Aip + f{t),ip{0) = ipo, (3.6) 

tffipit)^Tit)^o + J^Tit~s)f{s)ds. 

We call (j) G C'^([0, T], <¥) a weak solution to 

dtS{t)4>^A4> + f{t),4>{Q)^^o. 

iffv{t) = S{r-'^{t))(t>{r~'^{t)) is a weak solution to withf{t) = f{r-^{t))S{r-^{t)) 

and ipo ~ S{0)4>o- We call uj G C'^([0, T\,X^ a weak solution to 

S{t)dtLJ ^AU + f{t),Lu{0) = LUo, 

iff ip{t) = Lu{r^^{t)),ipo ^ ojo is a weak solution to h3.6]) . 

Remark 3.2. The given definition VS. 6\) is as the definition of weak solutions 
in Volume 5] or, equivalently, to the definition of mild solutions as given 
e.g. in f^. Since S is positive, the function r is invertible and due to the initial 
condition r(0) = 0. In the case S = 1 all definitions coincide. The motivation 
for the latter definitions is as follows: Provided we have sufficient regularity we 
can compute dt {S(t)(j){t)) = = A(f> + f{t) and, similarly, we compute 

S{t)dtuj{t) = dr^p{r(t)) = Auj + f{t). Similarly, a weak solution on [0,Tmax] x 

-dt4> = A4> + f{t), m^4>o 

is defined by a weak solution (pit) = ^(T'max — t) to i3.6\) with f{t) ~ f(T — t). 
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Next, we state semigroup properties for (|2.3p in the case S = 1. Denoting 
by C{X, y) the set of Hnear and bounded operators, under the assumptions 
(A1)-(A3), we have [TUlin] 

{J:ip) := at{x, n)ip e C {L^{Z x S'^),L^{Z x S'^)) , 
{Kip) = [ asix, il' ■ il)ip{x, n')di}' e C {L^{Z X S^),L^{Z x 5^)) 

and 

(Aifi) = -nv^^ip : D{A) L^{Z x 5^) 

is an unbounded operator on x 5^) with domain of definition D{A). We 

define 

T -.^ A~J: + K : D{A) C L^{Z x S^) L^{Z x S^) (3.7) 



and the general problem p.Sp reads 

d,Si' = T^ + q, i;^4^oeZ X X {0} and i/' = € r_ x [0, Emax]- (3.8) 



We summarize the properties of the defined operators in Proposition l3.31 see The- 
orem 2 in Volume 6, Chap. XXI and Definition 1 and Theorem 1 in Vol- 
ume 6, Chap. XXI for definition and properties of the trace of i/j G D{A). 
Due to Remark 3 in Volume 6, Chap. XXI. we have that ip G D{A) implies 
ip e L^(r_|ri • n{x)\dT^), where n{x) is the outer normal at a; G dZ. 

Proposition 3.3. Assume (A1)-(A3). Then, 

(1) A is the infinitesimal generator of a strongly continuous semigroup in 
L^iZ X S^). 

(2) T is the infinitesimal generator of a strongly continuous semigroup in 
L^iZ X S^) and D{T) = D{A). 

(3) For any %j.iQ G D{T) there exists a unique classical solution ip G C"'^ ([0, Cmax]; L^{Zx 
S^)) n C°([0, Emax]; -D(r)) of the Cauchy problem 

d,tp = TV, ^-(0) = V^o, - on r_ x [0, e„,ax] 

The solution is given by 'ijj{e) — exp{eT)'iJjo where cxp(eT) is the strongly 
continuous semigroup with generator T. 



For the non-homogenous problem (|3.8p and 5=1 the following result is 
classical [TJ[Tn] and in fact, it sufficics to verify the regularity properties of tp as 
defined below. 

Proposition 3.4. Assume (Al)-(A3), let q G L^([0, emax] x Z x 5^) and 
ipo G L'^{Z x S^). Then, for S = 1, there exists a unique weak solution ip G 
C°{[0,e^,^],L^{Z X S^)) to given by i' = exp(er)i^o + /q exp((e - s)T)q{s)ds. 

If additionally TpQ G D{T) and q G C^([0, Emax], -^^(-^ x 5^)), then ip is a clas- 
sical solution to (KW- Wehaveipe Ci([0, e^ax], i^(^xS'2))nCO([0, e^ax], £'(T)). 

If Q ^ ^ ^-''^d ^0 ^ 0, then we obain tJj > 0. 
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Proof. A proof of the second and third statement is given by Theorem 3 in 
Vohmie 6, Chap. XXI of [lU]- The first statement is classical and a proof is 
found Lemma 11.14, Theorem 11.16 in [H]. ■ 

Under the given assumptions on ipQ and q we do not necessarily obtain ?/' S 
C°([0,emax],£'(r))- However, if Vo S D(T) and q G L^{Z x S^), then V £ 
C^([0,emax],i^(2' X 5^)) n C°([0,emax],£'(T)). Additionally, wc have (Theo- 
rem 3, Chap. XXI, Volume 6, that if belongs to W2 then it is a 
pointwise a.e. solution to (j2.3p and it is unique in this space: For functions 
ip G W2 wc may apply Green's formula to obtain for any r G [0, emax] ■ 



(V(e),ifV'(e)-SV(e))rfe- 
llV-Wf L2(2x52) + (llV'(e)f L2(r_^) - ||V'(e)f L2(r_)) rfe > IIV'Mf l2(2x52) 

and hence \\iP{t)\\'^ l^zxs^) ^ 2 ||Al /J" ||'i/;(e)||\2(2xS2)^f and yields unique- 
ness due to Gronwall's lemma. The previous Proposition 13.41 also allows to 
define a control-to-state operator at first in the case 5 = 1 by 



Vo) =i' = cxp(eT)V'o + / exp((e - s)T)q{s)ds. (3.9) 
Jo 

Lemma 3.5. Assume (A1)-(A3) and let q G L^(Z x S*^ x [0, emax]) and ipQ = 
0. Then, the operator 'E, is a linear and bounded operator from L^{Z x x 
[0,emaJ)^i'(^x52 X [0,e^ax]). 

Proof. Note that for any strongly continuous semigroup exp(er) there exists 
a constant uj > and M > 1, such that || cxp(eT)|| < A/exp(ea;), see [50]. 5 
is linear and since C''([0, emax], -^^(^ x S"^)) C i^([0, emax] x Z x S*^) we have 

P(9)IIl2(zxs=x[o,6_.] ^ ^e'^'"";'' ■JI9IIl2(zxs=x[o,..„.J)- , ■ 
This yields existence of a minimizer for a class of cost functionals. We obtain 

the following result as extension to Theorem 3.1 in [TT] . 

Theorem 3.6. Assume (Al)-(A3). Let ai G L°^{Z) be positive a.e., let {[) G 

L'^(Z X X [0, eniax]) and q G L^^. Then, the problem 



/ 






/ 















subject to (3.10b) 



[ asi'dn' + q, (3.10c) 
iP ^OonT., iP{x,n,0) ^0 on Z X (3.10d) 

admits a unique minimizer q* G and corresponding weak solution ip* G 

C0([0,e^ax],i'(^x52)). 

Proof. For given i/j E L'^{Z x S"^) the operator 

$ : L2(Z X 52 ^ [0, e„,ax]) ^ ^'(^ X [0, e,nax]) (3.11) 

given by <&('!/') = ^g2 V' ~ ^/^df^ is an affine linear, bounded operator. Due to 
Lemma l3.5l the operator S is a linear and bounded operator from G L^{Z x 
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S"^ X [0, Emax]) L'^{Z X S'^ X [0, emax])- The subspace L^^ is a closed subset of 
a Hilbert space. Due to Theorem 2.16 m [57] the convex functional 

/ ^<^{Eiq)fdxde+ f ^{q-qfdxdnde 

attains its minimum. The latter being unique provided that a2 > 0. I 
Similarly, wc obtain an existence result for the equation with S ^ 1. 

Theorem 3.7. Assume (Al)-(A4). Let ai S L°°{Z) he positive a.e., let ip e 
L'^{Z X X [0, emax]) and q € L'^j^. Then, the problem 

min / — [ ~ i^] dxde + / ^^{q ~ qfdxd^de 

subject to 

dtS^) + QVxi' + cTtV = / (Jsipdn' + q, 

ip = Oon r_, tp{x, n,0) ^0 on Z X 

admits a unique minimizer q* S L^^ and corresponding weak solution ip* £ 
C0([0,e„,ax],i'(^x52)). 

Proof. Due to (A4), the solution -^r — ■5^,?'(0) = is a strictly mono- 
tone, smooth function with smooth inverse r~^. For any given q(x,fl,e) € L^^ 
introduce the operator ^ : L^^ by 

q{x,n,e) :=g(x,f7,r-i(e))5(r-i(e)) 

and we denote by (/) = E{q) for S defined in ()3.9p . Due to Proposition 13.41 6 G 
C°([0,emax],i^(^xS'2)) exists. Due to DcfinitionOV'(2:, f^, e) = (t)ix,^,r{e))/Sie) 
is a weak solution to equation (|2.3|1 with zero boundary conditions since r(0) = 0. 
We have -0 G C°([0, emax], -^^(^ x -S^)) and wc define the operator 

X:Ll,^ L\Z X X [0, e^ax]) : X{q){x, n, e) = ^i'l)) ^ (3.12) 

Due to Lemma 13.51 the continuity of S on the closed set [0, emax] and the 
C"'^— property of r. the operator X is a linear bounded operator. Then, the 
proof is exactly as in the Theorem 13.61 when replacing S by X. ■ 



4 First-order optimality conditions 

We define the operator 

T* := -A-J: + K : D{A*) C L^{Z x S^) L^{Z x S^) (4.13) 

for 

D{A*) := {A e L^{Z x S^) : nV^X £ L^{Z x S^), A|r+ = 0} 
and study for some r £ L'^{Z x S"^ x [0, emax]) the ec^uation 

- deX = T*X + r, X{x, emax) = 0, A = on r+. (4.14) 
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Lemma 4.1. Under the assumptios (A1)~(A3), there exists a linear and bounded 
operator E* : L^{Z x x [0, e„iax]) L^{Z x x [0, Cmax]) with X = E*{r) be- 
ing the weak solution to 1^. j^p for any r G L^{Z x x [0, emax])- Additionally, 
we have A G C°([0, emax], -^^(-^ x S*^)). T/ie operator E* is the adjoint operator 
on L'^{Z X S*^ X [0,emax]) to S provided that S, S* : L^(Z x 5^ x [0, emax]) — > 

W2nCO([0,emax],i?(AW)) CL2(ZX 52 X [0,emax]). 

Proof. Under the given assumptions and due to Theorem XX in Volume 6, 
Chap. XXI [ini, the operator T* is the infinitesimal generator of a strongly 
continuous semigroup with domain of definition D{A*). Hence, the equation 
= + / with zero initial data /^(O) = admits a weak solution given by 
IJi = /;exp((e - s)T*)f{s)ds £ C0([0, emax], x 5^)) for any / e L\Z x 

S"^ X [0, ema x]). H ence, A(e) = Ai(emax - e) G C°([0, Cmax], ^^(^ x S'^)) is a weak 
solution to (|4.14p . The solution A can be written as A(e) = — //"""L,, exp((emax — 
e — s)T*)f{s)ds. This defines a solution operator S* and as in Lemma [3.51 the 
operator S* is linear and bounded. Due to Theorem 3.3 in [TT] the operator T* 
is adjoint to T on x 5^ x [0, emax])- Given any z,w € L'^{Z x S'^ x [0, emax]), 
we denote by A = S*(z) G W2 and by ?A = S(u;) G W2. Note that for A, ■0 G W2 
we may apply Green's formula and hence obtain 



< W,5*(z) >L2(ZxS2x[0,e„ax]) = 



VXe),-f].V,A(e) + ^^^'^ 



de 



dt+ 

de = 



{Ki^{e) - S?A(e),A(e))^2(2xS2) 

(V'(emax), A(emax))L2(zxS2) ^ (^-(0) , A(0)) ^2 >< S^) + 

(^(e), A(e))p - {ijj{e),\{e))^_ de =< S(w),z >L^zxs^x[o,e,„^^]) ■ 



The previous result can be used to deduce the existence of a first-order 
optimality system in the case 5=1. This result will be extended to the more 
general case below. 

Theorem 4.2. Assume (Al)~(A3) and {j) G W2. Further, assume that the 
minimizer (q*,'ip* = S(g*)) of i3.10\) belongs to L^^ x W2. Then, the first-order 
necessary optimality conditions are 

q* = (q* - A* - a2{q* - q))+ a. e. e Z X X [0, emax], 

provided that A* = S*(ai Jg^ ip* - H^) and A* G W2 n C°([0, emax], -D(^*))- 

Proof. The proof is similar to the proof of Theorem 3.3 [TT]. Consider the case 
■0 = first. Then, the problem (j3.10p can be rewritten as 

miii|iy^*(2(9))|li2(2x[o,,„^j) + ll^/^('7-9)lii2(2xS2x[o,e„,,x]) =:.?(?) 

where $ is defined in (|3.1ip and S defined in p.9p . The adjoint operator on 
X [0,emax]) to $ : X 52 X [0,emax]) ^ ^'(^ X [0,emax]) Satisfies < 

W,$(v) >L2(Zx[0,e„„x]) = < **(w),W >Li{ZxS'^x[0,e^^„])= jzxS'^xyo^e^^^]^ ^ 

Due to Lemma 2.20, Lemma 2.21 and Theorem 2.22 of [37] and Theorem 3.3 
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in [TT], we obtain the gradient by j'{q) = {ai^{:i{q))) + a2{q - q). 

Under the given regularity assumptions and denoting by A* = S*(ai(?/'*) and 
"0* = S(g*) we rewrite the gradient as 

j{q) = S*(q!i / i^d^) + a2{q - q) = X + a2{q - q). 

In the case -0 we consider tp such that ^(yj) = ijj and as in Theorem 3.3 
we consider an optimal control problem for -0 = -0 — -0 subject to S(g — (p) = "0. 
Since S is linear and together with the result in the case -0 = we obtain the 
assertion. The strong form of the necessary conditions is given by system p.lOp 
in the case S = 1. ■ 
Proof of Theorem 12.21 . Existence of a minimizer is given by Theorem 13.71 
The minimizer q* G L^^^ is unique provided that a2 > and 0* is a weak 
solution in the sense of Definition 13.11 It is unique in the space W2 . 

Let r be the solution to r'(e) = 5^,''(0) = 0. Since S > and continuous, 
the function r is on [0, emax] and strictly monotone. Therefore, r is invertible 
and r' > 0. For q G L^{Z x S"^ x [0, Cmax]), and T defined by (031), let ip be the 
weak solution to (|2.3p . i.e., 

deS^p = T^ + q, n, 0) = 0, 1/' = on r_ . (4.15) 

Under the given assumptions and for any q G L^{Z x x [0, Emax]) there exists 
a weak solution ip G C°([0, emax], -^^(•Z' x 5*^)) C L'^{Z x x [0,emax]) unique 
in W2 due to Theorem 13.71 Then, 

cP{x,n,r{e))/S{e) ■.= ^{x,^,e), q{x,n, e) ^ q{x,^,r-\e))S{r-\e)) 

is a weak solution to 

= T(j3 + q, 0(x, n,o)^o, = on r_. (4.16) 

The solution operator is denoted by 2 as in p.9p . Similarly, for (0, q) being a 
weak solution to ()4.16p . there exists a unique ('0,'?) being a weak solution to 
(|4.15p . Hence, there exists a bijective mapping on C°([0, erna.^], [Z x S'^)) x L^^ 
from solutions(0, q) to (|4.16p to weak solution (0, q) to (|4.15p . The mapping pre- 
serves the regularity of the weak solutions, i.e., \i tp £ W2 fl C°([0, enmx], D{Aj) 
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then (/i e W2 n C°([0, emi,^],D{A)). Hence, we have for = r(e„ 



/ 



J 

J[o,. 



]xz 2 J 52 



— (q — (7)^ (a;, fl, e)dxdflde — 

[0,e„ax]xZxS2 2 



(4.17) 



xz 2 7^2 V S'(e) 

/* _ 2 

/ TT {q ^ q) {x, e)dxdVlde 



[0,Cmax]xZxS 



i^{x,n,r-'{T))dn 



012/ g(x,ri,T) _i \^ drdxdfl 



(4.18) 



where q are weak solutions to (|4.16p . Hence, if {ijj,q) of (|4.17p yields a min- 
imizcr {(t>,q) for (|4.18p . To obtain the necessary conditions, we may apply 
Theorem 112] to 

^° := Sir-\T))i:{x,n,r-\T),f = 5(r-i (T))g(a;, O, (r)) 



5(r-i(r)) 



since G L°°(Z x 5^ x [0,e,„ax]) and > a.e., G i^^^ ^ 5*2 x [0,e„iax]) 
and V'" e W2. We obtain on [0, Tr] x Z x S'^ : 

gO(a;, n, r) = i^q^x, n, r) - X°{x, il, r) - a2(|^^ - q{x, n, r^' (r))) ^ , 

A°(x,f},r) = S*(ai / fei^ _ ^(a;, 1], r^^ (r))df7), 
752 ^(r ^(t)) 



4>''ix,n,T)^E{q"{x,n,T)). 

Finally, we rewrite the optimality system in terms of and q. From the 
last equation we obtain that q{x,Q,,e) := q^{x,fl,r{e))/S{e) and ip{x,^,f-) ■= 
(j)°{x,fl,r{e))/S{e) is a weak solution to l|4.15p on [0, Cmax] x Z x 5*^. The sec- 
ond equation gives a weak solution to the following adjoint equation with zero 
terminal conditions and zero boundary conditions on r+ 

^ drX°{x,n,T) ^T*X"{x,n,T) +ai ( ^lj{x,n,r^^{T)) ~{l}{x,n,r-^{T))dn 



(4.19) 

Due to Remark [3.2l and Definition l3.1l we conclude from (|4.19p that A(a;, il, r^^(T)) := 
\^{x, r2, r) is a weak solution on [0, Emax] to 

-S{e)de\{x,n,€)=T*\{x,n,e)+ai I tP{x,n, e) - ij{x,n, e)dn. (4.20) 



12 



Further, A as the same regularity as A°. FinaUy, we reformulate the first equation 
of the optimality system. The equation is equivalent a.e. in [0, Tr] x Z x to 

Since S > and r bijective, we obtain using the functions tp and A for a.e. in 
[0, e„,ax] xZxS^ 

qix, e) - (q{x, fl, e) - A(x, Q, e) - ^2 {q{x, e) ~ q{x, fl, e))) ^ (4.21) 

The equations (|4.151 14.201 14.2ip comprise the first-order optimality system pro- 
vided that A, q and and therefore 0, A° and q'^ fulfill the given regularity 
assumptions. I 

5 Remarks and further discussion 

We offer the following remarks and notes on further discussion to the results of 
Theorem [ 



• The given results extend to minization problems of the type 

— [ [ '^fdxdVLde + — [ [ {q - qfdxdnde{b.2i) 

for given functions 9 S L^{Z x S'^ x (0, Emax)) by setting $ = Id. 

Under additional assumptions on the minimizer q* the assertions on ?/;* 
and A* can be obtained, e.g., if q* G W^'^{[0, e^nax], D{A)) we obtain that 
E{q*) G W2 n C°([0, e,nax], -D(^)) and similarly for A*. 

• Formally, we obtain the assertions of Theorem l2.2l as follows. Let (p{x, ft, r{e)) 
ip{x,fl,e) with r'(e) = and r(0) = 0. Then, the functional is 



[0,e,„ax]xZ 2 Jg-2 \ S{ 

f ^2 - 2 

+ / TT (q ^ q) {x, ^, e)dxdVtdt. 

"'[0,emax]xZxS2 2 

Using a coordinate transformation in e by r = r(e) and Tr = r(eniax) we 
obtain 



and its formal derivative in direction 6(f> as 



^« / /• o ^^ \ A^^a-^f]'^,-) dxdfl'dr 



ai 

JZxS^ 

Tr 



ai 

Jo Jzxs^ \Js- 



x,n,T) 



i}{x, vt, r-^ {T))dn ] 6(j){x, n', T)dxdn'dT. 



13 



• Instead of applying the transformation r'(e) = we could also study 
the properties of an energy-dependent advection operator. If we introduce 
V = i'fi^) with /(e) = l/5'(e) wc obtain the family of operator 

A{e) ■.= f{e)n-V,. 

Since / > this operator generates an evolution system G'(e, s) with 
domain of definition 

and by 

• In the beginning of Section[2]we assumed that at and as are independent of 
the energy level e. This is of course not a realistic assumption. Of course, 
formally, there will only be small changes to the first-order optimality 
system p.Sp in the case of energy dependent coefficients. One simply 
replaces at and as through their energy dependent counterparts. However, 
from an analytical point of view the semigroup theory presented here to 
solve the problem has to be extended in order to treat now evolution 
equations. Furthermore, the transformation e —>■ r(e) used in order to 
establish the main result cannot be used to simplify the problem. All 
these points will be discussed in future work. 
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